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O ■ 1 Problem 
o 

^ ' What is the shape of an asteroid of given uniform density and given total mass such that 

Qj • the force of gravity is maximal for one point on its surface? Compare the maximal gravity 

Q . with that of a sphere of the same mass and density. 

m 

^! 2 Solution 

Ph: ^ ^ ... . 

Jrj ■ This problem is from Yakov Kantor's Physics Quiz site, http:/ /star! tau.ac.il/QUlz7 



If the asteroid is spherical with radius a and mass M, the force of gravity on a test mass 
m is 



GMm 

everywhere on the surface, where G is Newton's constant of gravitation. 
r-| ■ Can some other shape of the asteroid result in a larger force? Let the test mass be at 

Ph| the origin, which is on the surface of the asteroid, and define the z axis to be along the 

direction of the desired maximal force on the test mass. It is "obvious" that this z axis is 
^ ■ an axis of symmetry of the asteroid: if the asteroid were not symmetric about the z axis, 

Q\ ■ the z component of the vector force Gmm'r/r'^ would be increased by moving material from 

^ . larger to smaller transverse distances from the axis, until the asteroid is axially symmetric. 

^ ! We seek the functional form x = x{z), where < 2 < ^o, that generates a surface of 

revolution about the z axis, the surface of the asteroid, that maximizes the (axial) force on 
• test mass m at the origin. The axial force due to a ring of extent dx dz that passes through 

■ point (x, z) is 

O ' 1 z 

' Gmp 2'Kx dx dz— , = , (2) 

^ . SO the total force is 

Ph. 

^= 2"^"^^ r " ^" (x^ +'.2)3/2 = r "^'i'- (x2(.) +.2)1/2 ) ' 



where we suppose that the asteroid has the same density as a sphere of radius a and mass 
M, i.e., p = SM/Aira^. Hence, the volume of the asteroid is constrained to have value 

V = 7T dz x\z) = . (4) 



This suggests use of the calculus of variations, although there is the complication that 
the value of the endpoint zq is unknown. 
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In principle, the shape x{z) might involve nonzero values at z = and 2:05 planar 
surfaces that bound the asteroid. However, we readily convince ourselves that such a geom- 
etry could not maximize the z component of the gravitational force; material in the plane 
z = contributes nothing to the z component of the force, while material at x > in the 
plane z = zq would be more useful if it were moved to a smaller z. Thus, we understand 
that x(0) = = x{zo) for the desired shape function. 



2.1 An Intuitive Solution 

However, we first consider the related question: what is the shape of a thin shell of matter 
such that all points on the shell contribute the same axial force on our test mass at the 
origin? For a shell described by x{z), the answer is already contained in the form of integral 
Namely, we desire that 

1 

const. = -2 5 (5) 



(a;2(2) +^2)3/2 



where we have evaluated the constant at the intercept zo of the shell with the z axis. On 
rearranging eq. (0), we have 

It is noteworthy that this form also intersects the z axis at the origin, so is a satisfactory 
candidate for the shape of the asteroid (whose surface surely must touch the test mass at 
the origin). Inserting this form in eq. (j3)) for the volume, we find that zq = = 1.7a. 
Using eq. © in eq. for the axial force, we find 

This argument does not prove that eq. (jH)) describes the shape of the maximal gravity 
asteroid, but it does show that the shape is not a sphere. The simplicity of the argument 
suggests that the result is in fact correct, but a proof of this is still desired. 



2.2 Solution via the Calculus of Variations 

We wish to find the form x{z) that maximizes integral F of eq. (0) subject to the constraint 
that integral V of eq. Q has the stated value, while permitting the endpoint zq to vary. To 
proceed via the calculus of variations, we consider the integral 

1= / f{z,x{z)) dz, (8) 

where ^ 

which combines the integrands of eqs. Q and (jH) using the Lagrange multiplier A. 
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For a fixed endpoint Zo, integral I is at an extremum provided / obeys the Euler-Lagrange 
equations 

df d df 

— = -. (10) 

dx dz dx' 

Since / does not depend on x' = dx/dz in this problem, we simply have 



^ ^ 9/ ^ xz 

dx (x^ + 2;2)3/2 



2\x. 



or 



2A. 



'11^ 



(12) 



(X2 + 22)3/2 

As in eq. (0), the constant 2A is clearly l/z^, since x{zq) = 0. Again, we obtain eq. © as 
the shape of the asteroid. And again, the constant zq {i.e., the multiplier A) must be set to 
v^a to satisfy the volume constraint (HJ, which completes the solution. 

The shape ^ of the maximal gravity asteroid is shown below. It is noteworthy that 
there is no cusp at either z = or zq. 



X 0.4 




2.3 A Mathematical Footnote 

We have not directly included in our calculus of variations the fact that the endpoint zq is 
free to move along the z axis, i.e., along the curve g{z) = 0. In, for example, sec. 4.3 of 
we read that the freedom to vary the endpoint in the maximization of integral / leads to the 
additional relation 

/ + (^'-^')^ = 0, (13) 

at the variable endpoint Zq, where g{z) describes the curve on which the endpoint zq can 
vary. Applying this to the present problem, we must have f{zo) = 0, which requirement is 
satisfied by the condition x{zo) = that has already been included in the analysis of sec. 2.2. 
However, it is pleasant to have a formal confirmation of our earlier "intuitive" justification 
of this condition. 
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2.4 Further Comment on the Shape the Asteroid 



The form ^ for the shape of a maximal gravity asteroid differs only slightly from that of a 
sphere, as shown in the figure on p. 3. Starting from a sphere, the force of gravity at, say, 
the "north" pole is increased by moving material from the "southern" hemisphere into the 
"northern" hemisphere. 

We verify that this is the right thing to do by a simple calculation. Namely, we compare 
the z component of the force of gravity at the north pole due to a point mass on the surface 
of a sphere of radius a at polar angle 6 with that due to a point mass at angle n — 6. The 
claim is that the point at 6 contributes more to the force than the point at tt — ^, so that it 
would be favorable to move some of the mass to 6 < 7t/2. 

The distance from the pole to the point at angle 6' is r = 2a sin ^/2. The z component of 
the force of gravity due to a point mass at angle 6 is proportional to 



Az a(l — cos6') 
7^ ~ 8a^sm^9/2 4:a^sm9/2 



^^W^^-S-^^TT.T^. (14) 



SO that the force due to a point at angle tt — 6* is related by 



The ratio, 

FJtt- 



=tan^/2, (16) 

is less than one for all 6' < tt/2. Thus, the force of gravity at the north pole would be 
increased by moving mass from angle vr — ^ to angle 6 on the surface of the asteroid, as 
claimed. 

While this argument gives a qualitative feel as to how a spherical asteroid should be re- 
shaped to become a maximal-gravity asteroid, it does not provide a quantitative prescription 
of how much mass should be moved. For that, the arguments of sees. 2.1 and 2.2 are better. 
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